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Abstract. We study the performance of quantum annealing for the simple p-body 
infinite-range ferromagnetic Ising model. In particular, we generalize the transverse 
antiferromagnetic interactions proposed by Seki and Nishimori as a quantum driver 
to many-body transverse interactions to understand if the two-body interactions 
are essential to allow the system to avoid troublesome first-order quantum phase 
transitions. We conclude that the general many-body interactions are effective to 
let the system evolve only through second-order transitions as long as a few minor 
conditions are satisfied. It is also discussed whether the overlap of the ground-state 
wave function of the new driver term with the target ground state is an essential factor 
for the success. 
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1. Introduction 

The task of finding the configuration that optimizes a given cost or energy function 
T-L({Si}) dependent on a large number iV of variables S±,...,Sn (often subjected to 
constraints) is a whole research field by itself, common to many fields in science. Finding 
the minimum energy or cost often becomes a hard task when the constraints in the 
system, or the interactions between variables, induce frustration because there is no way 
to find a minimum configuration that minimizes the problem locally. The frustration 
leads to a rugged landscape of many relative minima, and an exhaustive search for 
the absolute minimum is just not feasible for the interesting sizes (the dimension of the 
system often grows exponentially with N). As examples of these optimization problems, 
one can cite the traveling sales problem [Tj or the fc-SAT problem [2] in computer science, 
or finding the ground state of spin glass in physics [31 H] . 

Complexity in optimization problems is commonly classified as P if an algorithm 
is known to solve the problem in a time that grows polynomically with N. On the 
contrary, if it is not the case, and the time scales faster with N, these problems are 
labeled NP and considered as hard problems. Among all the NP problems, there is a 
subgroup named NP complete so that any possible NP problem can be reduced to one 
of them by means of a polynomial algorithm. Thus, if one algorithm were found that 
solved polynomically an NP complete problem, the whole family of problems would also 
become easy. The problems mentioned above belong all to the NP-complete classic 

Statistical mechanics, based on physical intuition, has contributed a lot in the 
development of new strategies for optimization problems: parallel tempering or replica 
exchange [6], and simulated annealing [7] are the two popular and widely used examples. 
In this last method, fluctuations are introduced in the problem through a ficticious 
temperature. This temperature favors the jump over barriers and thus encourages the 
system to visit other possible minima. The system is then simulated at a temperature 
T(t) that decreases slowly with time until it is finally switched off at the end of the 
simulation. We will refer to this simulated annealing as classical annealing (CA) in 
contrast to the quantum annealing (QA) [8j [9j [TOj [11], where fluctuations are induced 
also in the system but this time quantum ones. Quantum perturbations allow tunneling 
effects, and thus, if narrow enough, barriers can be crossed instead of surpassed. 

In the traditional QA formulation, a time- dependent Hamiltonian is introduced 

H(t) = s(t)H + [1 - s(t)] V, (1) 

where Ho is the target Hamiltonian (or the cost function that one wants to minimize) and 
V represents the quantum perturbations. In the field we are working in, the Hamiltonian 
Ho represents the magnetic interaction between spins. For the sake of simplicity, we 
will consider that H only depends on the z components of the Pauli matrix of, where 
i{= 1, . . . , N) labels the index of each spin in the system. As normally, we are interested 

| With the exception of the 2-SAT problem and the 2D Ising spin glass [5] that can be solved 
polynomically. 
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in finding the lowest energy spin configuration, i.e. the ground state. Now we introduce 
the quantum fluctuations through a spin driver term V. In principle, this term is 
arbitrary, as long as it does not commute with Hq. In addition, we impose that V 
has a single, trivial ground state. A typical example of a driver Hamiltonian is the 
transverse-field operator 

N 

VW = -J>f, (2) 

t=i 

where the of (i = 1, . . . , N) are the x components of the Pauli matrix. This perturbation 
is very intuitive, since it represents nothing but the interaction with a magnetic field 
along the x direction that induces quantum transitions between the eigenstates of of, 
whose modulus is tuned through the control parameter s(t). Initially, at t — 0, the 
control parameter s(t) starts at s(0) = 0, with H(0) = V, and increases monotonically 
with time until it reaches unity at time r and H{r) = Hq. Let us choose the simplest 
possible scheme where the control parameter grows linearly with time, i.e. s(t) = t/r. 
The evolution of the system, | $(£)), is determined by the Schrodinger equation, 

i^ |$( t )) = H (t) |*(t)> , 0<t<r. (3) 

The initial state |$(0)) is the ground state of the driver Hamiltonian V and is thus 
known. If the parameter s(t) is changed very slowly (r is very long), the state will 
be at every time very close to the instantaneous ground state. If it so, by tuning the 
parameters, one will move adiabatically from the initial ground state to the ground state 
of H . 

The adiabatic theorem states that the system stays close to the instantaneous 
ground state as long as r > A m f n where A min is the minimum energy gap from the 
ground state. Of course, in order for the above argument to be of general use, this 
A min cannot decrease with iV too fast. In fact, if the energy gap decays exponentially 
with the system size, as happens generally in first-order transitions, the running time 
will increase exponentially with N and the QA would not help to solve the problem 
efficiently 

This vanishing exponential gap present in many first-order transitions is sometimes 
considered to be one of the most important drawbacks of quantum annealing. Its 
presence was somehow shadowed for certain time by the preasymptotic behavior 
displayed in the small system sizes feasible in simulations (T2J E2 EH] • Indeed, in the last 
years, an increasing number of first-order transitions in the annealing parameters are 
being found [T5J EH HTJ [T8j EE] ■ It has thus been suggested that the presence of these 
quantum first-order transitions when tuning the transverse field is an intrinsic property 
of the systems with complicate free energy landscape, i.e. the hard problems, leading a 
pessimistic scenario for the QA algorithm [151 E3 CEO EIEl [19] . 

Recently, it was found that the ferromagnetic p-spin model, a model without 
disorder and with a simple free energy landscape, also suffers from this kind of first- 
order transition [18] . Due to its simplicity, this model constitutes a perfect benchmark 
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to study the QA performance. Indeed, it was recently shown [20j that, at least for finite 
values of p and p ^ 3, it is possible to avoid this first-order transition by appending an 
additional antiferromagnetic driver term and performing the annealing along a curve in a 
space of two annealing parameters instead of just one. This study changes the paradigm 
about first-order transitions in QA, since the failure of QA strategies observed up to 
now could be a failure of the standard formulation of QA with a transverse field, not a 
failure of the algorithm itself. 

Here we go deeper into this problem, studying a family of alternative driver terms, 
displaying different symmetries. We show analytically the existence of paths that cross 
only a second-order transition and thus the speed of QA is not exponentially damped. 
Indeed, in a second order transition the gap vanishes only polinomially with the number 
of particles, which must be compared with the exponential damping observed in the first 
order transition. The solution to the problem is not unique and we study the properties 
of these new driver terms, reaching the conclusion that the structure of the ground state 
of the additional Hamiltonians is not the main important feature that makes the whole 
algorithm success as argued in [2Tj . 

2. Problem 

Our starting point is the ferromagnetic p-spin model (p = 2, 3, 4 . . .) 




The ground state for this model, |$o)> corresponds to the state of all the spins aligned 
along the z direction. In order to avoid the degeneracy of the up and down configurations 
present in even powers of p, we consider here only the odd values of p and p > 3. In the 
limiting p — > oo case, this model is nothing but the Grover problem [181 122] • Although 
the Grover's quantum algorithm, whose reformulation in quantum annealing is given 
in [23], is considered a success of the quantum algorithm (provides a square-root gain 
with respect to the classical search [21]) it remains being a hard problem even with 
quantum algorithms. Now we consider the problem of finding this already known ground 
state |$o) of (jl]) with the QA algorithm using two driving terms. 
As usual, we consider the traditional transverse field operator, 

JV 

VW = -J>?, (5) 

i=l 

whose ground state, |$ TF ), is the one where all the N spins are pointing to the positive 
direction along the x axis. We next introduce a second Hamiltonian inspired in the 
antiferromagnetic interaction suggested in [20] . 




Many-body transverse interactions in the quantum annealing 



5 



that depends on a parameter k(> 1). When k = 2, we recover the antiferromagnetic 
interaction studied in [20]. The ground state for this Hamiltonian, namely depends 
on the value of the power k. When k is odd, the energy is minimum when all spins are 
aligned along the x axis but pointing to the negative direction. On the contrary, when 
k is even, the ground state corresponds to the state with total £) i=1 °t = if AMs odd, 
or YliLi a i = f° r N even. One of the goals of the present paper is to clarify whether 
the value k = 2 is essential to avoid the first-order transition. 

If we sum up (IT]), (EJ) and (jHJ), the new Hamiltonian of the problem reads as 

H(s, A) = s [XH + (1 - X)V k ] + (1 - s)Vtf- (7) 

Here there are two annealing parameters, s and A. These parameters will be tuned 
slowly during the annealing process so that, at the final time, r, s(r) = A(r) = 1 
and the target Hamiltonian (jlj) is thus recovered. In that way, one can explore the 
annealing process following infinitely different paths. It might resemble the idea of 
nondeterministic Turing machines, but one must always keep in mind that, even though 
many paths are possible, only one is chosen in each particular realization. 

The traditional QA is one of the infinite possible paths in (j7|). In fact, one can 
remove the influence of Vk, just by fixing A(i) = 1. Then, the annealing is performed 
by tuning s from to 1. If one looks at the configurations, at t — all spins should be 
aligned with the x axis, and at the end, with the z axis. In this case, we know that the 
system suffers from a quantum first-order phase transition between these two states. 
This transition ruins the efficiency of the algorithm as it becomes exponential [18J. 
The idea of introducing this two-parameter space (A, s) is precisely to try avoid this 
transition by following an alternative route. Seki and Nishimori succeeded in finding 
ingenious paths [20] with antiferromagnetic interactions, and here, we generalize that 
method to check how the value of k affects the conclusion. 

3. Analysis by a semi-classical approach 

The QA strategy will succeed if we are able to find a path in the space of parameters 
(A, s) that avoids crossing any first-order transition. With this aim, we compute in this 
section the phase diagram correspondent to the new Hamiltonian (EJ), as a function of 
the parameter k. The A" — > oo limit can be computed analytically using a semi-classical 
approximation (method to be explained below) or the Trotter- Suzuki decomposition 



formula [25] and the static approximation (see Appendix A ), leading to equivalent 
results. 

3.1. General Properties 

As a starting point, let us rewrite the Hamiltonian (J7J) in terms of the total spin variables 
(S a = | 5^i=i a i with a = x, y and z), 

2 ,Y , , / 2 x h 



H(s,X) = -s\N[-S z \ +a(l-\)Nl-S*) -2(1 -s)S*. (8) 
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This Hamiltonian commutes with the total squared spin, S 2 . Since the total spin is 
conserved and the initial state in the annealing process is the one with all spins aligned 
with the x axis, we are only interested in studying the maximum possible S value, i.e. 
S = N/2. 

Now, consider the normalized variables m a = S a /S, with a = x, y and z. The 

commutation relations for these variables are 

2 

[m x ,m y ] = l jj mZ i (9) 

and cyclic permutations. The normalized variable m a can take N + 1 values within the 
interval [—1, 1]. Thus, in the large N limit, these variables commute, and we can consider 
them as the components of a classical unit vector, i.e. m = (cos 9, sin 9 sin ip, sin 9 cos ip), 
being 9 the polar angle measured from the x axis, and ip the azimuthal one measured 
from the z axis. 

Considering the system now as classic, we can write the energy per spin as 

e = -s\(sm9 cos</?) p + s (1 - A) cos fc 9 - (l-s)cos0. (10) 

The equilibrium state will be determined by the minimum of e. Since p is odd, the 
minimum lies on the plane with ip = 0, which we call XZ + plane. The energy on this 
plane is labeled only by the polar angle 9 

e = -s A sin p 9 + s (1 - A) cos fe 9 - (1 - s) cos 9. (11) 

We search the #0 G [0, 7r] that minimizes ( flT]l l§l. The condition for the minimum is 

— = -psXsm p - 1 9 Q cos9 -ks{l-X) cos^" 1 9 sinfl +(l-s) sinfl = 0,(12) 
o9 

whose solutions are the angles 9q that satisfy either sin #0 = or 

ps\ sin p - 2 9 cos 9 + k s (1 - A) cos^ 1 9 - 1 + s = 0. (13) 

These two equations have more than one solution, and each one corresponds to a different 
phase. We will consider them as ferromagnetic if m z (= sin#o) > 0, and quantum 
paramagnetic if m z = 0. The most stable one at each point (A, s) will be the absolute 
minimum of e. 

We begin with the quantum paramagnetic solutions. The equation sin6*o = is 
satisfied for 9 = or it. The case 9 = corresponds to positive x magnetization, 
m x = 1. We name this phase QP + . Its energy is obtained by inserting this angle in 

eqp + (s,A) = s(l-A)-l + s. (14) 

The other paramagnetic solution, 9 = ir, corresponds to negative magnetization, 
m x = —1. We call this phase QP~. This phase is only stable for odd values of k 
and its energy is 

e QP -(s,A) = -s(l-A) + l-s. (15) 

§ Negative magnetizations along z axis have always higher free energy due to the change of sign in the 
sin p 8 term in (jTTJ) (remember that we only consider the p odd case in this work) . 
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This phase will not appear in the phase diagrams for k even, since its energy is always 
positive in the range of parameters < s, A < 1. 

We consider next the ferromagnetic solutions (8q > 0). The purely ferromagnetic 
solution sin# = 1 is only a valid solution on the line s = 1. Apart from this line, 
equation ( 1131) cannot be explicitly solved for any value of p, but it can be done in the 
p — > oo limit. We study below all the solutions for this limit and discuss their validity 
for p finite. 

3.2. Phase diagram for p — > oo 

In this limit, ( 1T31) has two possible ferromagnetic solutions. The parameter p appears 
in ffl3|) through psin p ~ 2 9 . We consider the two possible limits for the sine power, 
sin p ~ 2 6*o — >l (for the F phase) and (for the F' phase), always keeping 9 > 0. 
We begin the discussion with the F phase. With this aim, we assume 



In the p — > oo limit, this equation can only be satisfied if either the cosine vanishes, 
i.e. #0 = 7r/2 (but only on the line s = 1), or pcos9 tends to a constant. Let us 
investigate this second case. We consider cos#o = c/p, with c a ^-independent constant, 
and introduce it in (fTTl) . and taking the p — > oo limit, the equation reads 



whose solution is c = (1 — s)/s\. Thus, 

1 — s 

cos 6 = r^O, (19) 

sp A 

is a solution to (fl3|) . Still we need to check that this #0 agrees with the initial assumption 
(USD- Indeed, 



sin p # ->■ 1, 
and substitute it in (fl3l . 



(16) 



p s X cos 6»o + k s (1 - A) cos fe 1 9 - 1 + s = 0. 



(17) 



sAc-l + s = 



(18) 




We obtain the energy for this phase introducing (Tl9l in (fTT]) 




(20) 



On the other hand, the F' solution is obtained assuming the opposite limit, 



psin p 6 ->■ 0. 
Under this assumption, ([TBI) reduces to 

k s (1 - A) cos^ 1 6 - 1 + s = 0, 



(21) 



(22) 



whose solution is 



cos ft 



1 - s 



k-X 



(23) 



'0 = 



fcs(l-A) 
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Note that if k is odd, the negative solution for the cosine is also a valid solution. However, 
it has always higher energy than its positive counterpart, so we will not consider it for 
further discussions. 

The energy for the F' phase when p — > oo is 



e P /(s, A) | 



k - 1 



p— >oo 



1 - S 



l 

fc-i 



. M i-A)j ;i " s) - (24) 

Up to this point, we have obtained all the possible solutions to ( TTBl in the p — > oo 
limit: three (for even k) and four (for odd k) phases. We can use the energies to 
determine which phase is the most stable at each point (A, s). We show in figure 
[H several phase diagrams for k = 2, 3, 4 and 5. Let us analyze the nature of each 
transition. We begin with the transition line between the F' and QP + phases. This 
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line is obtained by solving e-pi — eQp = using the expressions (123)) and (JHJ). This 
equality is fulfilled on the line s = 1/[1 + k(l — A)]. On this line, m x = cos8 = 1 
in both phases, which corresponds to a second-order transition. On the other hand, 
the transition between the F and the QP + phases lies on the s = 1/2 line and, since 
magnetization is discontinuous, it is first-order. The second-order transition extends 
from (A, s) = (0, l/(k+ 1)) to (A, s) = ((k — l)/k, 1/2), the point where these two kinds 
of transitions cross. According to that, the higher k is, the broader the second-order 
line and the smaller the QP + region are. Furthermore, in the k — > oo limit, the QP + 
region completely disappears. 

Still there is a first-order transition between the F and F' phases, determined by 
the solution of ep> — e-p = using (1231) and (1201) . We solve this equation numerically and 
obtain the curve displayed in figure HJ On this line, the magnetizations are discontinuous 
but at the point (A, s) = (0,1) where they two become equal, m z = sm9 = 1. The 
transition is then first-order, but in the mentioned point, where it would be second-order. 

Up to this point, the discussion is common for even and odd values of k. However, in 
this latter case the QP~ phase also exists. Thereby, two additional transitions between F 
or F' phases and the QP~ phase appear. In both cases the x magnetization changes the 
sign on the transition, and then, they are first-order. The transition lines are obtained 
by solving the equations eq P - — ep = 0, leading to s = 1/(2(1 — A)), and e QP - — epi = 
which must be solved numerically. We display all the transition lines in figure HI 

According to these results, when we consider the p — > oo limit, there is only one 
single path that succeeds in avoiding first order transitions. This is the straight line 
that joins the initial point (A, s) = (0, 0) with the left upper corner, (0, 1), and the final 
state (1, 1). However, even though this path only crosses second order transitions, along 
this way there is no quantum annealing process, as can be seen by an insertion of these 
parameter values into the Hamiltonian ((7j), and thus this path is meaningless. 

4. Phase Diagram 

The phase diagram for finite p is different. Now, there appear regions where first-order 
transitions disappear, leaving more space for annealing trajectories. We display the 
corresponding diagrams in figures El 0, SI and [5] for k = 2, 3, 4 and 5, respectively. 
Again, the shape of the phase diagram strongly depends on whether k is even or odd. 
In the former, there are only three phases and in the latter the extra QP~ phase appears. 
Besides, the higher k is, the longer is the second-order transition line. 

The picture of the ferromagnetic phase for finite p is rather complicated. When 
one solves numerically ( ITBl and looks at the 6 Q > solutions, the situation is the 
following: in a wide region, one finds two possible alternative solutions that look very 
much alike to the F and F' phases discussed for the p — » oo limit. However, near the 
left and upper corner in the phase diagram, there is one single ferromagnetic solution 
which is neither F nor F' but something intermediate. In fact, for k even, one can 
find paths through which the magnetization evolves continuously from the F' to the F 
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Figure 2. Phase diagram for k = 2. This is the same phase diagram as in reference 
|20) . The transition between the F' and QP phases is of second order, and the F-QP 
and F-F' transitions are of first order. 
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Figure 6. Magnetization obtained with the semi-classical approach as a function of 
s for A = 0.1 and for k = 2 and 3. The dashed lines correspond to the analytical 
predictions for the QP ± , F <(26]) and F' ([32]) solutions. 



magnetizations without crossing any transition on the way, see figure [6j However, when 
p is high and k is odd, transitions between the F and F' phases cannot be avoided, see 
figures |6] and [71 

All this effect can be understood quantitatively coming back to the discussion of the 
p — > oo ferromagnetic solutions. Each of the phases were derived using the assumptions 
( TT6|) for the F phase, and (l2Tj) for the F' phase. Now we discuss the validity of these 
approximations for p finite. 

We begin with the F phase. This phase was obtained by introducing ffl6|) in f fl3|) . 
Since this equality is not strictly true, we introduce it as an approximation sin p ~ 2 6*o ~ 1, 
thus obtaining a new approximate equation 

p s A cos 6 + k s (1 - A) cos^ 1 9 - 1 + s w 0. (25) 

If we assume cos^o 1, the solution is 

cos6 ^ sipxl'kii-^y (26) 

for k = 2, and 

cos # ~i^i (27) 

Sp A 

for k > 2. That means, that the F solution found for the p — > oo limit also appears for 
finite p whereas cos#o <C 1, or 

^— f « 1- (28) 



Many-body transverse interactions in the quantum annealing 13 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

s s 



Figure 7. Magnetization obtained with the semi-classical approach as a function of s 
for A = 0.3 for k = 2 and 3. The dashed lines correspond to the analytical predictions 
for the QP ± , F {26} and F' (|32|) solutions. 



In particular, the smaller this quotient f[28|) is, the better approximation the F solution 
is. We can obtain the energies for finite p by introducing this solution in ( TTTT) . For 

k = 2, 



e k F = 2 (s,\); 
+s (1 - A) 



-sA 



1 



s[pA + 2(l-A)] 



s[pA + 2(l-A)] 

A- 

- (1 - S. 



1 - s 



s[pA + 2(l-A)] 



and for k > 2 



ep(s, A) w -s A 
1 - 



1 - 



1 



+s(l-A) 



A 



sp A 
fl- 



1 - s 
sp A 



i 

fc-i 



(29) 



(30) 
(31) 



Next we study the F' solution. We consider the following approximation 
psin p - 2 e w 0. 
As before, if this is a good approximation, 

cose °^[k7^X}_ 

is one solution to (|T3|) . This solution is equal to the one obtained for p — > oo, (|23l) . In 
other words, at this order of approximation, the solution is exact at this limit. 



(32) 
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We briefly discuss the range of validity of this F' solution ( 13 2 p for p finite. The 
approximation ( l3~lj) is valid for small values of 6*0. With this idea we expand separately 
the two terms in ( Tl~3l) around 8q = 0, the first term being 



sin p 2 #o cos 6o = p s X 8q 2 

and the second term 

ks(l - X)cos k ~ l e - 1 + s 



l-^o 2 + O(0 o 4 ) 



fcs(l-A) 



i - k —±el + o(e { 



-l + s. (33) 



The dependency on # in the first term becomes irrelevant when p > 3, thus recovering 
the F' solution (|32|) . When p = 3, the lowest power of 8 appears in the first term, 
leading to a different ferromagnetic solution, but not the F'. Clearly, the higher p (and 
the smaller # ) is, the better is approximation fl3T|) . 

In general, for intermediate values of s and A, the higher p is, the more exact the 
two ferromagnetic solutions, F and F', are. Then, since both approximations represent 
opposite cases in the value of m x (or m z ), a new first-order transition between both 
phases will appear on the line when their two free energies become equal. However, for 
low values of p, or alternatively for s — > 1 or 0, there will only be one ferromagnetic 
solution, somewhere in between these two F and F' phases. This idea is well ilustrated in 
figures [6] and [7J where both the numerical solution to (fl"3l and the analytical predictions 
(1261) and (1231) are displayed. 

This has straightforward consequences on the performance of the quantum 
annealing algorithm: the higher p is, the narrower will be the region where annealing 
paths can avoid a first-order transition. In the limit of p — > oo, as was discussed before, 
there will be only one possible path, but not effective as quantum annealing. 

Concerning the transitions between the QP and ferromagnetic phases, we can 
distinguish two kinds of transitions. First of all, the transitions between the F and 
QP phases will be first order, since the F phase is characterized by a high value of m z 
whereas the paramagnetic solution has m z = 0. On the other hand, there is another 
transition between the F' and QP phases that lies on the line where their two free 
energies become equal, i.e. s — 1/[1 + k(l — A)]. On this line, m x = 1 (m z = 0) for the 
two phases. Furthermore, the F' solution is exact for m x = 1. Since the magnetizations 
are continuous on this line, the transition between F' and QP is of second order. Besides, 
it can be checked that there is a wide range of this line where ep> < e-p. Thus, this phase 
is the stable one in the ferromagnetic phase. This second-order transition does not 
hamper the QA performance and gives us a way to avoid the F-QP phase transition 
that appeared when using the traditional QA approach. It is important to point out 
that this second-order transition appears for any value of k. 

In Appendix A , we describe a different, quantum-mechanical method to derive the 
same results. 
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5. Energy gap 



As discussed in Introduction, the efficiency of the QA algorithm is closely related to 
the behavior of the gap between the ground and first excited states. As usual, this gap 
can be computed by direct diagonalization of the problem Hamiltonian (|7j). Indeed, 
since the total spin S is conserved during the evolution, the dimension of the problem is 
N+l. That means that the diagonalization matrixes grow polynomially with the system 
size instead of exponentially as for generic quantum problems. However, still computer 
resources limit this computation to moderate sizes although such computations are 
useful for some purposes (2TJJ, [18]. Here, we adopt an alternatively approach, this 
gap can be computed in the thermodynamic limit iV — > oo by the method described 
in [26]. The main idea is to extend the semi-classical scheme for the ground state 
by the consideration of quantum fluctuations around the classical ground state. It is 
important to point out that this method can only be applied in the case of finite gaps in 
the thermodynamic limit, as it is the case away from the transition points themselves. 
In case of exponentially small ones, other methods such as instantonic or WKB methods 
should be used [HI EI]. 

It is most convenient to rotate the system by an angle 9q around the y axis in order 
to bring the x axis parallel to the semi-classical magnetization, i.e. 



(34) 
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We rewrite the Hamiltonian (jSj) in terms of these new variables S a , obtaining 



H(s,X) = -s\N 
+ s(l — X)N 



2 

N 



cos 9 S x + sin 9q S z 



2 

N 



2(1 



sin 9 S x + cos 9q S z 
sin 9 S x + cos 9 S. 



(35) 



Now, we add quantum fluctuations to the system by means of the Holstein 
Primakoff transformation 
N 



S z 



— a) a 



S, 



(N 



a) a) 1 ' 2 a 



~Sl 



where a is a boson annihilation operator that satisfies [a, a'] = 
fluctuations are small relative to the classical state, i.e. for iV ^> 
simpler expression 

/N 



S x 



(a + a 



(36) 

1. When quantum 

(aJa), we can use a 

(37) 



We introduce these transformations into the Hamiltonian (|35[) and expand the three 
different terms in powers of 1/N. Thanks to the previous rotation, the coefficient in 
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l/y/N vanishes. We keep terms up to l/N and group together all the coefficients with 
the same power of N. The result is 

if( 7> 5) = jVe + 7 + 7 [(a f ) 2 + a 2 ] + 5a) a. (38) 
The term for iV 1 is nothing but the ground energy obtained before in (|lip . 

e = -sX sin p 9 + s(l — A) cos fc 0„ - (1 - s) cos 0„. (39) 

The coefficients S and 7 are given as 

5= — sX [p(p - 1) sin p ~ 2 fl cos 2 O - 2p sin p 0„] 
+ s (1 - A) [Jfe(Jfe - 1) sin 2 6 cos k - 2 6 - 2 k cos k O ] + 2(1 - s) cos O ,(4O) 

and 

7 EE _ SA ^~ 1 ) sin P-2 £ Q CQS 2 £ Q + s (1 _ A ) fc (fc-1) gin 2 ^ CQS k-2 6o{41) 

We need to diagonalize this Hamiltonian in order to compute the first excited state by 
the Bogoliubov transformation 

a = cosh — b + sinh — b\ 0) = cosh — b^ + sinh — 6, (42) 
2 2 2 2 v y 

where b is a new bosonic annihilation operator satisfying [b, & '] = 1. Using this 

transformation, we can eliminate the coefficient of [(b^) 2 + 6 2 ] by choosing the angle 

as 

tanh 6 = — J- = e. 


With this choice, the Hamiltonian can be written as 

5 



H(j, 5) = Ne + 7 + - (Vl - e 2 - l) + A (43) 



with 



A = c^VT^e 2 . (44) 

The Hamiltonian is diagonal in tfb. The energy gap in the iV — > 00 limit between the 
ground and first excited states is A. 

Using the values 9 previously obtained solving (fi~3j) . we can compute the energy 
gap for our system. We show the data for p — 11 and A = 0.1 and 0.3 for different values 
of k in figure El As was suggested in the magnetization data in the previous section for 
A = 0.1 (figures 121 to EJ) , no first-order transition F-F' is observed through the energy 
gap. The gap vanishes continuously on the second-order transition line but present no 
further jumps later, but the ones related to the F-QP - that always take place in the 
odd-fc cases. On the contrary, when A = 0.3, the jumps in the gap appear for all the fc's 
at the place where we observed the F-F' transition before. 

In the thermodynamic limit, the gap vanishes at a single point of first-order 
transition and remains finite away from this point. The single point of vanishing gap is 
hard to see by the present method, which results in an apparent simple jump in the gap 
at a first-order transition as seen in figure El 
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6. Overlap of the ground-state wave functions 

It has been suggested in [21] that the reason for the antiferromagnetic interaction, the 
k = 2 case in (jH]) introduced by Seki and Nishimori in |20j, to work better than the 
transverse field interaction only is related to the large overlap between the ground states 
of the Hamiltonians V^=2 and H . In this section, we will discuss the properties of these 
different states, concluding that, even thought the overlap is important, it is not the 
decisive factor that makes the strategy to succeed. 

The ground state of Vtf is the one where all the spins are aligned along the x axis, 
|0tf) = ®fLi - If we denote the ground state of H , as \<p ) = |t)^, the overlap 
between |0tf) and |0o) decreases exponentially with N as 2~ N , as can easily be seen 
from the elementary relation = (|t)^ + \\)*) /a/2- 

The overlap computation becomes a little more complicated for the ground state of 
Vfc. The ground state for this term depends on the value of k. Indeed, if k is odd, the 
ground state is the one where all the spins are aligned along the x axis, but towards the 
negative direction, i.e. ®f =i |v)f • Then, the overlap with |0 O ) for the k odd case will be 
exponentially suppressed as 2~ N as in the case of Vtf- Thus, the argument in [21] does 
not apply directly since we can avoid first-order transitions even in this case of k odd, 
in spite of the very small overlap of the ground state for H and V^. 

The ground state for the k even case needs some care to be analyzed. We compute 



it in Appendix B We show there that that the overlap is indeed higher for k even. The 
antiferromagnetic interations is a particular case, k = 2. In fact, the overlap displays 
an algebraic decay as the system size increases, i.e. ~ 1/y/N. 

We conclude that the overlap is not the main ingredient that makes the present 
method to succeed. 
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7. Conclusions 



We have analyzed the reason for the failure of the traditional annealing with a transverse- 
field term in the infinite-range ferromagnetic p-spin model. We have shown that it is 
possible to find annealing trajectories that avoid the crossing of first-order transitions 
thanks to the introduction of a second driver term in the problem, which may be due 
to the multiple spin flips in the z-basis caused by the second term as was the case 
in [27]. This additional term favors the appearance of a second-order transition that 
does not hamper the annealing performance. A whole family of possible candidates has 
been studied and we conclude that the solution to the problem presented by Seki and 
Nishimori [20] is a special case of a more general additional quantum term. The main 
properties of these additional terms have also been discussed with the conclusion that 
the properties of the ground states of the diverse terms in the Hamiltonian are not a 
decisive factor to make the quantum annealing fail or succeed. 
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Appendix A. Analysis with the Suzuki- Trotter formula 



We investigate the properties of H(s,X), the phase diagram in particular, using the 
decomposition formula [25] and the static approximation. This approach, although 
quantum, leads to the same results as the semi-classical method described in section [3j 
The method here is analogous to the one explained in detail in [20], [18] , but we leave 
the power k as a free parameter in all the calculus. The purpose of this appendix is to 
confirm consistency between the method of the main text and references [201 IT8] . 
The starting point is the partition function, 

Z = Tre-^ (s ' A) . (A.l) 

We use the decomposition formula to express it as 

M 



lim Z M 



lim Tr 

Af-s>oo 



■^sAffog-^ [s (1-A)V A ff + (1- S )Vtf] 



lim > 

M-s-oo 



(K}|^exp 



x exp 



(3s (1 - X)N 
M 



(3s\N 
M 



N ^ 



0~i 



+ 



(3(1 -s) 
M 



N 



M 



IK}), 

(A.2) 

where Yl{a z } re f ers to the summation over all the 2 N possible spin configurations in the 
z basis, and |{cr z }) = ®f = i \of). 
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We introduce M closure relations, each one labeled by a(— 1, . . . , M), 

i(a)= IK(«)))(K(«)}|x E (K(«)}l , (A.3) 

just before the ath exponential operator involving in ( 1A.2I) . The trace over the 
product of quantum operators is thus reduced to the product of numbers that commute 
and can be reordered, 

M 

n e e «p 



Af 



x exp 



/3s\N 



M \ N 



N 



l N \ i" 



i=l 



08(1 -\)N ly^*. 



M 



a) 



+ 



(3(1 



i=l 



M 



N 

E 



of (a) 



N 



xll«( a ) l»f («)> (of(a) W(a + 1)>, 



(A.4) 



i=l 



where \af(M+ 1)) = 

We write the product in terms of the total x and z magnetizations in each copy of 
the system, i.e. m x (a) = jj J2iLi a i( a ) an d m z {a) = j?J2iLi a i( a )i using the integral 
definition of the delta distribution 

TV 



V i=i 



dm 5 mfa) 



1 N \ 

i=l / 



(A.5) 



After a few simplifications, we introduce the static approximation to remove the a 
dependence of the magnetizations. Under this approximation, we can compute the 
M — > 00 limit using again the decomposition formula. The partition function ( 1A.ll) 
then reduces to 



dm 2 dm x exp [-N/3 f(/3, s, A; m 2 , m x )} , (A.6) 

where /(/3, s, A; m 2 , m x ) is the pseudo free-energy defined as follows: 

f{(3, s, A; m 2 , m x ) = (p - 1) 5 A(m 2 ) p - (Jfe - 1) s (1 - A)(m x ) fc 



1 





-^log <{ 2cosh/3J[psA(m z )p- 1 ] 2 + [1 - s - s (1 - A) k (m 1 )^ 1 ] 2 J> (A.7) 



Again, one can apply the saddle-point method, obtaining two self-consistent equations 
for the two magnetizations, 

p s A (m z ) p ~ 1 



m 



^\ps\{m z )p- 1 } 2 + [1 - s - s (1 - A) (m x ') fc_1 ] 2 

x tanhPyJ\ps A (m 2 )^ 1 ] 2 + [1 - s - s(l - A) (mf 4 ] 2 , 
1 - s - s (1 - A) jfe (mf 4 

^/[psA (m 2 )^ 1 ] 2 + [1 - s - s (1 - X)k(m x ) k - 1 ] 2 
x tanh/3^/[ps A (m 2 )?- 1 ] 2 + [1 - s - s (1 - A) k (m x ) k - 1 ] 2 . 



(A.8) 



(A.9) 
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In this work we are only interested in the purely quantum transitions, 
not in the thermodynamical ones. For this reason, and with the sake of 
simplification, we remove the dependence of physical quantities on (3 from now on 
by considering the low-temperature limit, j3 — » oo. In this limit, if [ps A (m z ) p ~ 1 ] 2 + 
[l — s + s (1 — A) k (m x ) k ~ 1 ] 0, the hyperbolic tangent in (1A.8jl and (1A.9I) tends to 
unity, and thus the self consistent equations simplify 

m z = psXjmr- 1 (A 1Q) 

-,x\k-l] 2 



\j\ps\ (m z y- l f + [1 - s - s (1 - A) k (m 1 



m* = l-.-.(l-A)t(m-)>-' (A n) 

y[ps\ {m z yp- 1 } 2 + [1 - s - s (1 - A) k {m^y- 1 } 2 

The magnetization lies on the unit radius circumference, i.e. (m x ) 2 + (m z ) 2 = 1. This 
result agrees with the approach in section El where the magnetization was a unit vector 
constrained to the XZ plane. The pseudo free energy ( 1A.7j) becomes 



™* = , r fc u , (A.i3) 



f{fi, s, A; m z , m x ) = (p - 1) s X(m z ) p - (k - 1) s (1 - AXm*)* 

- ^/[j9 s A (m^P- 1 ] 2 + [1 - s - s (1 - A) k (m x ) k - 1 } 2 . (A.12) 

Equations ( 1A.10I) and ( lA.lip have ferromagnetic (F) solutions with m z > and quantum 
paramagnetic (QP) ones satisfying m z = and m x ^ 0. Let us begin with the latter 
case. 

Appendix A. 1. Paramagnetic solutions 

Substituting m z = in ( 1A.11I) . we get 

ks (1 — \)(m x ^ k ' i 
\l-s-ks(l- \){m x ) 

which leads to m x = ±1. The solution m x = — 1 is obtained if the numerator in (1A.13I) 
is negative, that is, if 1 — s — k s (1 — A)(— < 0, which, in the range of parameters 
< s < 1 and < A < 1 considered, can only be satisfied if k is odd and in the region 
1/[1 + k(l — A)] < s < 1. This phase is precisely the QM~ phase discussed in the text. 
Its free energy is 

/qp-(s,A) = 1-2s + sA, (A.14) 

which coincides with equation ( ITS]) . 

The other quantum paramagnetic solution with m x = +1 (the QP + phase) can be 
satisfied only if the numerator is positive, i.e. if 1 — s — k s (1 — A) > 0, which can be 
fulfilled for any value of k as long as s lies in the region <s < 1/[1 + k(l — A)]. The 
free energy of this phase is 

/ QP+ ( S ,A) = -1 + 2 S - S A, (A.15) 

and is also equal to (Q3J). 
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There is still one additional paramagnetic solution. In order to obtain it, we need 
to come back to the discussion about the /3 — > oo limit. The hyperbolic tangent in flA.8j) 
and ( 1A.9I) could tend to a finite value in the — > oo limit, as long as the term in the 
square root vanishes. Mathematically, 



when 



lim tanh (3\ [ps X (m 2 ^ 1 ! 

/3— >oo 



m z — > 0, m x — > 



+ [1 - s - s (1 - A) k (m 1 )"] 2 = tanhc, (A.16) 



(A.17) 



_Jfes(l-A). 

In order to find a non-trivial solution, it is also necessary in this limit that m z tends 



(A.18) 



to zero faster than the bracketed term of m x in flA.9j) . i.e. 

ps X(m z ) p ~ 1 

- — > 

l-s-ks(l- \){m x ) k - 1 

Under these assumptions, flA.8j) and flA.9j) imply m z = and m x = tanhc, where 
tanhc = [(1 — s)/ks (1 — A)]* 771 , in order to be consistent with the limit (IA.17p . This 
correspondence determines the region in the space where this phase can appear. In fact, 
as any hyperbolic tangent, | tanhc| < 1, which is true only if 1/[1 + k (1 — A)] < s < 1. 
Besides, the condition (IA.18P forces p > 30 

Since the magnetization in the z direction vanishes, we call this phase QP2. The 
free energy is obtained with (IA.7p . 



/qp2(s, A) 



k - 1 

k 



1 - s 



:i-s) 



(A.19) 



ks(l — A) 

This last phase was not predicted by the semi classical approach. However, we will 
see below that it is irrelevant to the problem, since the F' phase has always a smaller 
value of the free energy. 



Appendix A. 2. Ferromagnetic solutions 

We next consider the possible solutions with m z > 00 As before, the ferromagnetic 
solutions cannot be computed explicitly for a given value of p but for certain limiting 
cases. 

I In the fc-odd case, the limit 



0, m x 



fes(l-A) 



also makes the square root in (|A.16[) vanish, but it leads to a positive free energy in (|A. 19|i . and thus 
it is not relevant. 

Indeed, using (m x ) 2 + (m z ) 2 = tanh 2 c = [(1 — s)/ks (1 — A)] ^ and computing the limit (|A.17|) 
when m x — > tanhc, one can check that it vanishes only as long as p > 3. 
+ No negative value for m z can satisfy (IA.10[) for odd values of p. 
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The solution m z = 1 (and m x = 0) is exact only on the line s = 1. However, we 
can see that an approximate solution m z « 1 and m x ~ is valid in a wider space of 
parameters. Indeed, the solution 



m 



1 - s 
sp A 



and m z 



1 



spX 



(A.20) 



fulfills (lA.lOj) and (lA.llj) when (1 — s)/ps\ — > 0. This is the F phase we obtained 
before in equation ( 1301) . The free energy is obtained plugging these values into equation 
(IA.12D . For the p -)■ oo limit, 

f F (s,\)\ P ^co = -s\. (A.21) 

We consider an alternative solution for < m z < 1. With this aim, we rewrite 
(lA.lOj) in the following way 

[(mf - 1] [ps A(m*)* -1 ] 2 +{m z [l - s - s(l - X)k(m x ) k ~ 1 ] } 2 = 0.(A.22) 

In the p — > oo limit, p(m z ) p ~ 1 — > 0, and 



m 



1 - s 



_ks(l-X)_ 



m 



1-s 



ks(l - A) 



2 

fc-1 



(A.23) 



is an exact solution to (IA.22J) . and similarly of ( lA.llj) . as long as (1 — s)/ks(l— \)<s< 1, 
or l/[l + k (1 — A)] < s < lo This is precisely the F' phase discussed in sectional Again, 
we compute the free energy by plugging the solution flA.23j) in flA.12j) and taking the 
p — > oo limit 



Ms, A) | 



k-1 



p— >oo 



(A.24) 



_ks(l - A) 

which is exactly equal to the one obtained for the QP2 phase ( 1A.19I) . 

The solution (IA.23j) is also a good approximate solution for p finite (but p > 3) 
when (m z ) p — > 0. The free energy for this phase is 

2 TP 



Ms, A) 



-sA 



1 - 



1 



sk(l- A) 



k-1 

k 



1 - s 



_ks(l-X)_ 



(1-s), (A 



which, for finite p, is always smaller than /qp2- According to this observation, except 
for the p — > oo limit, the F' phase is always stabler than the QP2 phase. 

We have therefore reproduced the results of section [3] by a completely different 
method. The present method is nevertheless better suited for generalizations to more 
complicate problems where the target Hamiltonian Ho cannot be expressed in terms of 
simple total spins. 

* Again, the negative solution for m x is also a valid solution in the odd k case but has a higher free 
energy than (|A.23|) due to the change of sign in the (m x ) k term in (|A.12[) . 



Many-body transverse interactions in the quantum annealing 



23 



Appendix B. Ground state of V k and its overlap with the ground state of H 



In this Appendix, we derive the properties of the ground state of Vk for k even. Let us 



first consider the case with N even. The ground state of H , \<f> ) = 
expressed as 

I0o) = ®£i(lt>? + U>f)/V2 

= ( m • • • m + m m ■ ■ ■ m-, m 

+ --- + \i) X 1 \i) X 2 ---\i) X N 



1 |t) - , can be 



(B.l) 

This last expression has 2 N terms, in which the partial sum of terms with a half of the 
sites having |t)^ and the other half \l) x is nothing but the ground state of in the 
S = N/2 sector up to a normalization, 



10*) = a |t)? |t>2 • • • tt)w ID 



N/2 \+/ N/2+1 



ID 



N 



+ 



+ UK 11)2 ' ' ' \i) X N/2 It) 



N/2 I I /N/2+1 



lt&). 



(B.2) 



It is easy to check from the number of terms in the above equation that the normalization 
N \ 

„ T . =1. We thus have 
K N/2j 



condition is a 2 



(0o|0fc) = 



( N 



2 N/2 y N / 2 




For large N, 



log | (00 1 (p k ) | 2 = log 



-N 



(W. 



ilogiV + logJ| 



(B.3) 



(B.4) 



which means that the overlap decreases only polynomically with N &s ~ N 1 / 2 . 



or 



The case of odd can be analyzed similarly with 

' N 



N 
N/2 



replaced by 



N 

(N + l)/2 
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